Persistent random walk on a one-dimensional lattice with random asymmetric 

transmittances 
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We study the persistent random walk of photons on a one-dimensional lattice of random asymmet- 
ric transmittances. Each site is characterized by its intensity transmittance t (t' ^ t) for photons 
moving to the right (left) direction. Transmittances at different sites are assumed independent, 
distributed according to a given probability density J-{t,t'). We use the effective medium approxi- 
mation and identify two classes of T{t, t') which lead to the normal diffusion of photons. Monte Carlo 
simulations confirm our predictions. We mention that the metamaterial introduced by Fedetov et al. 
[Nano Letters 7, f 996 (2007)] can be used to realize a lattice of random asymmetric transmittances. 
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I. INTRODUCTION 

Random walks in random environments is a field of 
continuous research [l|, 0, H, H| • Hopping conduction of 
classical particles or excitations [Bj], transport in porous 
and fractured rocks |6j , and diffusive transport of light in 
disordered media 0, are a few examples. 

Random walks with correlated displacements figure in 
a multitude of different problems. Among the correlated 
walks, the persistent random walk introduced by Furth 
Q and Taylor Q , is possibly the simplest one to incorpo- 
rate a form of momentum in addition to random motion. 
In its basic realization on a one-dimensional lattice, a per- 
sistent random walker possesses constant probabilities for 
either taking a step in the same direction as the immedi- 
ately preceding one or for reversing its motion [3, H, 0] • 
Generalized persistent random walk models are utilized 
in the description of polyme rs [loll, chemotaxis [TTj ] . gen- 
eral transport mechanisms 12L 13]. Landauer diffusion 
coefficient for a one-dimensional solid [3], etc. 

Recently, the persistent random walk model is used to 
study the role of liquid films for diffusive transport of 
light in foams [TH, [l6| ■ Diffusing- wave spectroscopy ex- 
periments have confirmed the photon diffusion in foams 
(l7l |. A relatively dry foam consists of cells separated by 
thin liquid films [181 ]. Cells in a foam are much larger 
than the wavelength of light, thus one can employ ray 
optics and follow a light beam or photon as it is trans- 
mitted through the liquid films with a probability t called 
the intensity transmittance. This naturally leads to a 
persistent random walk of the photons. In the ordered 
honeycomb (Kelvin) foams, the one-dimensional persis- 
tent walk arises when the photons move perpendicular 
to a cell edge (face). Thin-film transmittance depends 
on the film thickness. Films are not expected to have the 
same thickness. These observations motivated us to con- 
sider persistent random walk on a one-dimensional lattice 
of random transmittances UM- We assumed that trans- 
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mittances at different sites are independent random vari- 
ables, distributed according to a given probability den- 
sity f(t). Assuming that < 1/t >= J Q f(t)/t dt is finite, 
we validated the classical persistent random walk with 
an effective transmittance t e ffj where l/t e // =< l/t>. 
We also investigated the transport on a line with infinite 
< 1/t >. We showed that if f(t) -> /(0) as t -> 0, the 
mean square-displacement after n steps is proportional 
to n/bx{n). If f(t) ~ f a t-° (0 < a < 1) as t -> 0, 
we found that the mean square-displacement is propor- 
tional to n( 2 ~ 2a ^( 2 ~ a K Quite interesting, we found that 
anomalous diffusion of persistent walkers and hopping 
particles on a site-disordered lattice 0, are similar. 
To observe photon subdiffusion experimentally, we sug- 
gested a dielectric film stack for realization of a distribu- 
tion f(t) [3. 

In the realm of diffusion on one-dimensional lattices 
with random hopping rates Wjj> from site j' to site j, 
the asymmetric hopping model with w -j -j+i ^ Mj+i.i 
has gained much attention [H EH, El El EI El SI : At 
variance with the symmetric case, the asymmetric model 
can display anomalous diffusion behavior without broad 
distribution of hopping rates. The asymmetric hopping 
model has been used to discuss hopping conductivity in 
pre sence of an external electric field, molecular motors 
[271 ] . evolution of a domain wall in a one-dimensional 
random field Ising model d , helix-coil transition of het- 
eropolymers pjl [2q. etc. These points suggest us to in- 
vestigate persistent random walk of photons on a one- 
dimensional lattice with random asymmetric transmit- 
tances. For a given dielectric stack, the transmittance 
for incidence on the right side, is equal to that for in- 
cidence on the left side [29]. However, optical elements 
with different transmission in the forward and backward 
directions, and even optical diodes which allow unidi- 
rectional propa gat ion, are realized [3f| . For example, 
Fedetov et al. [3l| showed that asymmetric transmis- 
sion through a planar metal nanostructure consisting of 
twisted elements can be observed in the optical part of 
the spectrum. For a normally incident circularly polar- 
ized light of wavelength 630 nm, this metamaterial is 
1.3 times more transparent from one side than from the 
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other. There is a good reason to believe that the exper- 
imental observation of photons' persistent random walk 
is not out of reach: Barthelemy, Bertolotti, and Wiersma 
have recently verified Levy flight of photons in their syn- 
thesized Levy glass j32j . 

Apart from interest in the optics of random media, our 
work has been motivated by the Lorentz gas model intro- 
duced to describe the diffusion of conduction electrons in 
metals [HI, [34[ . One-dimensional persistent random walk 
and stochastic Lorentz gas are intimately related (35|. 
Stochastic Lorentz model consists of fixed scatterers on 
a lattice and one moving light particle. The light particle 
runs at velocity c or — c, and when collides with a scat- 
terer it is transmitted with a site-dependent probability 
t or reflected with a probability 1 — t. In other words, the 
light particle performs a persistent random walk. Thus 
here we are investigating a variant of the one-dimensional 
Lorentz gas, where each scatterer is characterized by a 
random asymmetric transmission coefficient. 

In this paper we consider the persistent random walk 
of photons on a one-dimensional lattice of random asym- 
metric transmittances. Each site is characterized by 
its intensity transmittance t (t') for photons moving to 
the right (left) direction. Transmittances at different 
sites are assumed independent, distributed according to 
a given probability density J-(t, t'). We generalize a vari- 
ant of the effective medium approximation introduced by 
Sahimi, Hughes, Scriven, and Davis [ijj HH to identify 
two classes of T(t, t') which lead to the normal diffusion 
of photons: (i) (1/t) is finite and (t'/t) is less than 1. (ii) 
(1/f) is finite and (t/t') is less than 1. Here ( ) denotes 
averaging with respect to the distribution T(t, t'). 

Our paper is organized as follows. In Sec. |H] we in- 
troduce the model. In Sec. IIIII we present the effective 
medium approach to the problem. The numerical treat- 
ment and its results are reported in Sec. |TV] Sec. |V] is 
devoted to a discussion of our results. 



II. MODEL 

We consider a one-dimensional lattice random walk in 
which steps are permitted to the nearest neighbor sites 
only. We normalize the length and duration of a step to 1. 
Apparently, on a one-dimensional lattice the walker can 
move either to the right (+) or to the left (— ) direction. 
Each site j is characterized by forward and backward 
transmittances tjj+\ and t'- respectively: On arriv- 
ing a site j, a walker moving in the right (left) direction 
takes a step in the same direction with the probability 
tjj+i (tj Here we assume asymmetric transmit- 

tances, i.e. tjj+i ^ t'jj-v 

We assume that (i) transmittances t and t' at each site 
are random variables. In general, these random variables 
are not independent, (ii) transmittances at two differ- 
ent sites are independent, (iii) transmittances at all sites 
are distributed according to a given normalized proba- 
bility density T{t, t'). Apparently J Q J Q !F(t,t')dtdt' — 



1. The probability density functions of t and t' are 
/+(*) = Jo f{t,t')dt' and f~(t') = T{t,t')dt, respec- 
tively. The joint probability distribution can be written 
as T(t,t') = f + (t)f~(t') when random variables t and 
t' are independent. For any function h(t,t'), we define 
< h(t, t') >= Jq 1 Jo h(t, t')F(t, t')dtdt'. 

We denote by P + (n 7 j) (P _ (n,j)) the probability that 
the walker after its nth step arrives at site j with positive 
(negative) momentum. A set of two master equations can 
be established to couple the probabilities at step n + 1 to 
the probabilities at step n: 

P+(n+l,j) = tj^P+inJ -l)+r' j _ 1J p-(n,j -1), 
P-(u+1J) = r j+lij P + (n,j + l)+t' j+1J p-(n,j + l), 

(1) 

where Tjj-i = 1 — tjj+i and r'j . +1 = 1 — tj j_ l denote 
forward and backward reflectances at site j, respectively. 

We are mainly interested in the probability that the 
photon arrives at position j at step n, i.e. P(n,j) = 
P + (n,j) + P~(n,j), from which we extract the first and 
second moments after n steps as the characteristic fea- 
tures of a random walk: 

«J»» = (E^'i))' 
3 

({f))n = (£j 2 P(n,j)). (2) 

3 

Here the first bracket represents an ensemble average over 
all random transmittances, and the second bracket signi- 
fies an average with respect to the distribution P(n,j). 

Assuming a constant forward transmittance t and a 
backward transmittance t' at each site, translational in- 
variance of the medium can be invoked to deduce the 
exact solution of P(n,j) in the framework of charac- 
teristic functions Q. Furthermore, the mean square- 
displacement of photons after n — ► oo steps can be ob- 
tained as (j 2 ) n — (j)n — 2Dn, where the the diffusion 
constant D is 

n _ 2{\-t)(\-t'){t + t') 

U ~ (2-t-t'f ' [6) 

In the limit t = t' one obtains D = t/(2 — 2t), a known 
result in the realm of the the classical persistent random 
walk. 

The disorder not only may affect the value of diffu- 
sion constant as compared to the ordered system, but 
also may lead to the subdiffusive or superdiffusive be- 
havior. In our model, even a few sites with small trans- 
mittances may drastically hinder the photon transport: 
In the extreme limit where at two different sites j and j', 
transmittances = tjj-i — tj'J'+i = tj'j'-i = 0, 

photons either do not visit the segment between j and j , 
or are caged in this segment. In the following section, we 
determine which distributions of transmittances J-(t, t') 
lead to the normal diffusion of photons. 
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III. EFFECTIVE MEDIUM APPROXIMATION 



where 



Many of the approaches to the transport in disor- 
dered media have the disadvantage of being restricted 
to one-dimensional problems. Here we adopt the effec- 
tive medium approximation (EMA) which is applicable 
to two- and three-dimensional media. We generalize a 
variant of effective medium approximation introduced by 
Sahimi, Hughes, Scriven, and Davis [19l. l36||. 

First we simplify the set of coupled linear difference 
equations ([T]) using the method of the z-transform 0, l37| 
explained in Appendix [Al 



P+(z,j) P+(n = 0,j) 



P-(z,j) P-(n = 0,j) 



= tj-ijP+izJ-l) 
+r' j _ 1>j p-(z,j-l), 

= r j+hj P + (zJ + l) 

+t' j+1J p-(zJ + l).(4) 



We assume the initial conditions P + (n = 0, j) = P~ (n — 
Q,j) = Sjfi/2. To facilitate solution of Eq. (H|) we in- 
troduce probabilities P e ± (z, j), and a reference lattice or 
average medium with all forward transmittances equal to 
t e (z) and all backward transmittances equal to t' e (z), so 
that 



P+(z,j) P+(n = Q,j) 



P c 7(z,j) p-(n = 0,j) 



t e (z)P+(z,j -I) 
+r' e (z)p-(z,j-l), 
r e (z)P+(z,j + l) 
+t' e (z)P e -(z,j + l). (5) 



Here r e (z) = 1 — t e {z) and r' e (z) = 1 — t' e {z) denote 
effective reflectances. 

EMA determines t e (z) and t' e (z) in a self-consistent 
manner, in which the role of distribution J-(t, t') is man- 
ifest. This is done by taking a cluster of random trans- 
mittances from the original distribution, and embedding 
it into the effective medium. We then require that av- 
erage of site occupation probabilities of the decorated 
medium duplicate P^(z 7 j) of the effective medium. We 
will sketch the method in the following. 

Subtracting Eqs. ^ and (O, we obtain 



UQ + (z,jj 

z{Q-(z,j), 
't 



( JQ + (*J - i)\ _ T+( JQ + (zj + 1; 



3-1,3 r -- 




^i)-T-(z) 











-T+(z) 



P+(z,j-l) 
P-(z,j-l) 

P+(z,j + l) 
P~(z,j + 1) 



,(6) 



Q+(z,j)\ [P + (z,j)\ (P+(z,j) 
Q-(z,j) - \P-(z,j) {p-(zJ)J> 



T-(z) 



T+(z) 



t e (z) r' e (z) 






To solve Eq. ([B]), we introduce the Green function 

Gil G\2 



(7) 



G(zJ) 



G21 G22 



which satisfies the equation 

-G(z,j)-T-(z)G(z,j -l)-T+(z)G(z,j + 1) = S jfi I. 

z 

(8) 

Here I is the identity matrix. Multiplying both sides of 
the above equation by e lje and then summing over all 
the sites, the Fourier transform of the Green function i.e. 
G(z, 0) = e^GHz, j), can De obtained as 



G(z,0) 



A(z,B) \ r e (z) 



t' e (z)e-* 6 r' e (zy 



,-10 1 



t e (z)e l 



(9) 



where A(z,0) = f - z[t e (z)e ie + t' e (z)e- ie ] + z 2 [t e (z) + 
t' e (z)-l]. 

For the present, we consider only the simplest approx- 
imation, and embed one random transmittance at site I 
of the effective medium. Then solution of Eq. ([6]) is 



Q + (^T 
Q-{z,j) 



J " G(z,0)S{z,O)e- ie(j - l) 



where 



[t lil+1 -t e (z)]e* e [t' e (z)-t'u-i 



[t e {z)-t w ]e-^ [tU-i-t' e (z)}e-* >- {11) 



S(z,d) = 

Self-consistency equation is < P ± (z,l) >= P e ± (z,Z), 



< [I - J 2 J G(z, 0)S(z, 9)^-} - 1 >= I. (12) 

The above matrix equation leads to these conditions: 

pl rl T(t,t r )dtdt' 

J l-(t-t e (z))U(z)-(t'-t' e (z))U'(z)~ ' 

1 r 1 tF{t,t')dtdt 

Jo 1 - (t t e (z))U(z) (t> t' e (z)W'(z) te[Zh 
1 rl t'Fit^dtdt 1 _ , 

l-(t-t e (z))U(z)-(t'-t' e (z))U'(z) * e(z) ' 

(13) 



Jo 
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where 
U(z) = 

U\z) = 
V(z) = 



V(z)(l + z 2 (-t e (z) + t' e (z) -!))-! 
2t e (z) 

V(z)(l + z 2 (t e (z)-t' e (z)~l))-l 
W e (z) 

1 

+ z\t e (z) + f e (z) - l)) 2 - 4zH e (z)t' e (z) 

(14) 



It turns out that one of the self-consistency conditions 
(fT"3|) can be trivially satisfied. Consistency equations de- 
termine t e (z) and t' e (z), in which the role of distribution 
T(t, t') is manifest. For symmetric transmittances where 
T(t,t') — f(t)S(t — t'), our consistency conditions (TT3|) 
indeed duplicate that of Ref. [IB] 

The translational invariance of the effective medium 
can be invoked to access the z-transform of the the first 
and second moments of the photon distribution: 



n=0 



[j)nZ n 



t e (z)-t' e (z) 



n=0 



(i-z) 2 i-z[t e (z) + t' e (z)-iy 

2z 2 {t e (z)~t' e (z)Y 

(i-zf (i-z[t e (z)+t> e (z)-i]y 

Z l + z[t e (z)+t' e (z)-l] 



(i-z) 2 i-z[t e (z)+t e (z)-iy 



(15) 



We are interested in the long time behavior, thus Taube- 
rian theorems suggest to analyze Eqs. (|13[) and (|15[) in 
the limit z — » 1. 

To find which distributions of transmittances ^(t, t') 
lead to the normal diffusion of photons, we assume that 
t e (z) and t' e (z) have no singularity in the limit z — * 1. 
We find two distinct classes: (i) If if! ft) < 1 then 



t'M = 



i 



(I/*)' 
(t'/t) 
(i/t) ' 



The first class of admissible distribution !F(t, f') is such 
that (1/t) is finite and (t'/t) is less than 1. (ii) If (t/t') < 
1 then 



*e(*) = 

t' e {z) = 



Ml 
(i/ty 

i 



(17) 



The second class of admissible distribution T(t, f') is such 
that (1/t') is finite and (t/f) is less than 1. The above 
effective transmittances do not depend on z, thus one 
can directly use Eq. ([3]) to access the diffusion constant 
of photons. 

The Cauchy-Schwarz inequality states that for two 
random variables a and (3, (a[3) 2 ^ (a 2 )(0 2 ). With 
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FIG. 1: The diffusion constant D as a function of trans- 
mittances t and t' of an ordered medium. Theoretical and 
Monte Carlo simulation results are denoted, respectively, by 
lines and points. 



a = yft/V and /3 = ^JFJi, we find that 1 ^ (t/t')(f ft). 
This clearly shows that (t/f) and (t'/t) are not simulta- 
neously less than 1, thus two mentioned classes are quite 
distinct. EM A does not predict any result when both 
(t/t 1 ) and (t'/t) are greater than 1. 



IV. NUMERICAL SIMULATIONS 

The predictions of EMA can be inspected by numer- 
ical simulations. The computer program produces 50 
media, whose transmittances are distributed according 
to a given distribution T{t, t'). We deliberately focus 
on cases where both (t/t') and (t'/t) are not simulta- 
neously greater than 1. For each medium, the program 
takes 10 4 photons at the initial position j = and gen- 
erates the trajectory of each photon following a standard 
Monte Carlo procedure. The statistics of the photon 
cloud is evaluated at times n G [10000, 12000, ...,68000]. 
(16) ((j 2 ))n — ((j))n is computed for each snapshot at time n, 
and then fitted to 2Dn+0 by the method of linear regres- 
sion. An offset O takes into account the initial ballistic 
regime. We compare our numerical diffusion constant 
with the analytical one based on Eqs. ©, (fTB"|) and (JT7J). 

First we assume that transmittances of all sites 
are equal. For t e [0.3,0.5,0.7] and t' e 
[0.0, 0.1, 0.2, 0.9, 1.0], our numerical and analytical 
predictions for the diffusion constant are compared in 

Fig.m 

Next we consider F{t,t') = f + (t)f~(t') such that 
f + (t) is a uniform distribution for t\ < t < t2, and 
f-(t') = 5(t' - t' 3 ). We choose t x = 0.2, 0.2 < t 2 < 0.6, 
and tg G [0.1,0.7,0.9]. Our numerical and analytical 
predictions are compared in Fig. Uta). We also consid- 
ered the case where f + (t) is a uniform distribution for 
t\ < t < t2, and f~(t') is a uniform distribution for 
4 < f' < t' 4 . We choose t x = 0.1, 0.2 < f 2 < 0.6, f 3 = 0.6 
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FIG. 2: (a) The diffusion constant D as a function of t% for 
different values of t' 3 . T(t,t') = f + (t)f~(t'), such that f+(t) 
is a uniform distribution for 0.2 < t < ti and f~(t') = 5(t' — 
t' 3 ). (b) D as a function of £2 for different values of t' 4 . f + (i) is 
a uniform distribution for 0.1 < t < ti, and /~ (t') is a uniform 
distribution for 0.6 < t' < t'4. Theoretical and simulation 
results are denoted, respectively, by line and points. 



FIG. 3: The diffusion constant D as a function of a. J-(t, t') = 
f + (t)f-(t') such that (a) f+(t) = (1 - a)t _a for < * < 1, 
and f~(t') is a uniform distribution for 0.7 < t' < 0.9. (b) 
f + (t) = (1 - and /~(t') = 5(t' - Theoretical 

and simulation results are denoted, respectively, by line and 
points. 



and t' 4 6 [0.7,0.9]. Our results are shown in Fig. HJ^b). 

We also present two other examples. We consider 
F(t,t') = f + (t)f-(t') such that f+(t) = (1 - a)t- a 
for < t < 1, and f~(t') is a uniform distribution for 
0.7 < t' < 0.9. Fig. |3fa) depicts D as a function of 
a. Our results for the case f + (t) = (1 — a)t~ a and 
f-(t') = 5(t' - t[) are illustrated in Fig.^b). 

Figures [T][3] vividly show that the effective medium ap- 
proach to the diffusion constant D is quite successful. 



V. DISCUSSIONS 

In the present paper, we address the persistent ran- 
dom walk of photons on a one-dimensional lattice of ran- 
dom asymmetric transmittances. Clearly, a photon steps 
back by each reflection, fntuitively, one expects the abun- 
dance of large reflectances (e.g. two different sites with 
tjj+i = t'jj-i = tj'j'+x = t'j'.j'-i = 0) to drastically 
decrease excursion of the photons. As percolation prop- 
erties [l|, Q , this feature is induced by the dimensionality 
of the lattice. We focus on determining distributions of 
transmittances F{t, t') which lead to the normal diffusion 
of photons. The probability distribution P ± (n, j) as an 
exact solution of the master equation (fj| is quite hard to 



obtain. However the relatively simple but approximate 
effective medium approach reveals intriguing aspects of 
the system. In two cases, the transport of photons is 
diffusive: (i) is finite and (t'/t) is less than 1. (ii) 

(1/f) is finite and (t/t') is less than 1. Monte Carlo sim- 
ulations confirm our predictions. EMA does not predict 
any result when both (t/t 1 ) and (t'/t) are greater than 1. 

It would be instructive to compare our problem with 
the transport on a one-dimensional lattice with random 
asymmetric hopping rates. Independent steps are the 
base of the hopping transport, while correlated steps are 
the essence of the persistent random walk. Hopping con- 
duction is described by the master equation 

8P(t i) 

g T 1 = w jJ+1 P(r,j + 1) + Wj j-iP(T,j - 1) 

-Oi+l,i +Wj-lj)P(T,j), (18) 

where P(r, j) is the probability for the particle to be on 
site j at continuous time r, and Wjj> denotes the prob- 
ability of jumping from site j' to site j per unit time. In 
the asymmetric hopping model Wjj+x ^ Wj-\-ij. First 
we note that EMA does not predict any result when both 

(""M+iM+ij) = (w^/wj) and JtUj+ijMj+i) = 
(w^,/Wi-) are greater than 1 [23|, [26|]. Making use of 
a periodization of the medium, Derrida obtained exact 
expressions for the velocity and diffusion constant [2l| . 
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In the case (\og(w^/w^)) < 0, he found (i) The velocity 
V vanishes if (w^/w^) ^ 1. (ii) For (w^/w^) < 1 < 
((w^/w^) 2 ) the velocity is finite but the diffusion coef- 
ficient is infinite, (iii) For ((w^/w^) 2 ) < 1 both V and 
D are finite. All these results are easy to transpose when 
(log(u;<_ /«;_►)) > 0. 

As already mentioned in Sec. U the metamaterial in- 
troduced by Fedetov et al. [3l| can be used to realize a 
lattice of random asymmetric transmittances. We sug- 
gest a simple arrangement where a fraction e of the lat- 
tice sites are randomly occupied by the metamaterial, 
and the rest of lattice is occupied by half transparent di- 
electric slabs. One can measure the diffusion constant of 
photons as a function of e to test our predictions. For 
the proposed distribution of transmittances T(t, t') = 
e6(t - 0AS)6(t' - 0.57) + (1 - e)S(t - 0.5)6(1/ - 0.5), we 
find 



D = 



2(1 - 0.2456e) 
(2 -0.2456e) 2 ' 



EMA does not predict any result when both (t/t') 
and (t'/t) are greater than 1. It would be interest- 
ing to investigate the anomalous diffusion of photons 
and self-averaging quantities of the system following 
Refs. [Ill, [zj [23| . Our studies can also be extended to 
higher dimensional lattices. 
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APPENDIX A: 2-TRANSFORM 

The z-transform F(z) of a function F(n) of a discrete 
variable n = 0, 1, 2, ... is defined by 



F(z) = J2F(n)z n . 



(Al) 



One then derives the z-transform of F(n + 1) simply as 
F(z)/z — F(n = 0)/z. Note the similarities of this rule 
with the Laplace transform of the time derivative of a 
continuous function @, [3?J • 

Under specified conditions the singular behavior of 
F(z) can be used to determine the asymptotic behavior of 
F(n) for large n (Tauberian theorems) Q. For example, 
the identity T(l - a)(l-z) Q - 1 = T(n-a+l)z n /n\ 

shows that 



F(z) ~ m ~ Q) - F(n) = r( -" + 1) , (A2) 
where T(a) — J Q e~ t t a ~ 1 dt. Particularly, 



F(z) 
F(z) 



1 



(l-z) 2 
1 



F(n) = n + 1, 
1 



(l-z) 



, * F(n) = -(n z + 3n + 2). (A3) 
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